H 4 and other hydrogen-oxygen compounds at giant-planet core pressures 
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Water and hydrogen at high pressure make up a substantial fraction of the interiors of giant 
planets. Using ab initio random structure search methods we investigate the ground-state crystal 
structures of water, hydrogen, and hydrogen-oxygen compounds. We find that, at pressures beyond 
14 Mbar, excess hydrogen is incorporated into the ice phase to form a novel structure with H4O 
stoichiometry. We also predict two new ground state structures, P2\/m and I4/mmm, for post- 
C2/m water ice. 



As one of the most abundant substances in the solar 
system, water ice at high pressure [TJ [2] is of fundamen- 
tal importance in planetary science [3]. Over the last 
two years, a significant effort has been devoted to finding 
the ground-state structures of ice at the multimegabar 
pressures corresponding to the cores of gas giant planets 
(15-20 Mbar for Saturn and 40-60 Mbar for Jupiter). Mil- 
itzer and Wilson |4j demonstrated that the Pbcm phase 
[2] becomes unstable above 7.6 Mbar and predicted a 
new structure of Pbca symmetry. Wang et al. [5] showed 
that, at 8.1 Mbar, water ice assumes another structure 
of IA2d symmetry instead of transforming into a Cmcm 
structure as proposed in Ref. [4]. Refs. [6]-[8] instead pro- 
posed a structure with Pmc2i symmetry to appear in the 
same pressure regime, but it has a higher enthalpy than 
the lA2d structure. Around 12 Mbar, a transition to a 
structure with P2\ symmetry was consistently predicted 
in Refs. [5-8]. Even higher pressure studies predict a 
P2i/c structure at 20 Mbar [7], and a metallic C2/m 
structure at 60 Mbar [3 [8] . In nature, however, high- 
pressure water phases are rarely found in isolation, and, 
in gas giants, an icy core may be surrounded by a vast 
reservoir of hydrogen-rich fluid or exist in a mixture with 
the other planetary ices, such as ammonia and methane. 
Recent works have emphasized the fact that counterin- 
tuitive stoichiometrics can occur in post-perovskite ma- 
terials at extreme pressures [9]. This raises the question 
of whether H2O is indeed the ground-state stoichiometry 
for water at high pressure in hydrogen-rich environments, 
or whether pressure effects are likely to result in the in- 
corporation of hydrogen into the ice lattice to form novel 
ice-like phases with non-H^O stoichiometry. 

In this Letter, we apply ab initio random structure 
searching (AIRSS) methods [10] to determine the ground- 
state structure of hydrogen-oxygen compounds at ex- 
treme pressure to explore whether the H2O stoichiom- 
etry of water is maintained at high pressure. The AIRSS 
method relies on the generation of a large number of 
random geometries whose structures are then optimized 
using density functional theory. In the AIRSS process, 
randomly generated unit cells are filled with randomly 



positioned atoms, and the structures are geometrically 
relaxed to the targeted pressure. Structures with com- 
petitive enthalpy are picked out and re-evaluated with 
more accurate thermodynamic calculations, from which 
the most stable structure can be determined. Although 
the AIRSS is not guaranteed to find the most stable 
phase, it has achieved remarkable success in discovering 
structures across a wide range of materials pTtfl3] . 

In our calculations, structural optimization is per- 
formed using the Vienna ab initio simulation package 
(VASP) [14]. Projector augmented wave (PAW) pseu- 
dopotentials [15] and the Perdew-Burke-Ernzerhof [T6] 
exchange-correlation functional are used. The pseudopo- 
tential cutoff radii equal to 0.8 and 1.1 Bohr for H and 
O. For efficiency of the AIRSS, a cutoff energy of 900 eV 
is chosen for initial relaxation, which is then increased 
to 1700 eV for precise calculations on favorable struc- 
tures. Likewise, a relatively sparse Monkhorst-Pack grid 
[TTJ [18] (6x6x6) is used initially, which is later replaced 
by a 20 x 20 x 20 grid. Tests at the largest energy cutoff 
and grid size show that the total system energies are con- 
verged to the order of 1 meV, while the energy difference 
is converged to 0.1 meV. 

As an initial step, we recomputed enthalpies of re- 
cently reported [IHZ] H2O structures as a function of 
pressure, assembling all proposed structures for com- 
parison. Figure [I] shows the transition sequence up 
to the highest reported pressure (70 Mbar) as follows: 
X Pbcm Pbca Il2d P2 1 P2 1 /c -+ C2/ra, 
consistent with previous reports. Note that the Pmc2\ 
structure predicted in Refs. [6H5] is less stable than the 
l\2d structure reported by Wang et al. [5] at all pres- 
sures; as the Ii2d structure has eight formula units (f.u.) 
per unit cell, it was less likely to be generated by random 
search algorithms. 

With quickly increasing number of known massive exo- 
planets, it is worthwhile to explore stable H2O structures 
at even higher pressure. We use AIRSS to generate be- 
tween 1000 and 3500 H 2 unit cells with 1-4 f.u. at 
pressures of 100, 150, 300, 400, and 500 Mbar. Addition- 
ally, we generate between 750 and 1150 structures with 
5-8 f.u. at pressures of 100 and 200 Mbar. 

Figure [l] shows two new, stable phases with P2\/m and 



2 




-0 
-0 

|-0 

^ 0. 

Q. 

CD 0. 

_c 

c -0. 
LU 

-0. 

-0. 
-0. 



.00 
.05 F 
.10 
.15- 
.20-- 



15 




W^lceX[1] 
-OPbcm[2] 
Pbca[4] 
■-V-M2d[5] 
— ft^P21[5-8] 
Pmc21[6-8] 



00 f) — © — ®^© — ^-©-©©©-g-g ©©V©-®. 

15*-<C ' C2/m 

1 b >^ P21 [5-8] 1 VVV V P21/m ^ 
30-^VP21/c[7] ^V^.** 
45 :^P>C2/m[6,8] ^WVA-2U I 



-0-P21/m 
60 ^-V-l4/mmm 



20 40 60 80 




120 180 240 300 360 420 480 
Pressure (MBar) 



FIG. 1. (color online) Enthalpy versus pressure of different water ice structures. The reference following the name of each phase 
is the corresponding earliest report. In the lower panel, structures are provided for C2/m, P2\/m and 74/rarara, where large 
red and small blue spheres represent O and H atoms, respectively. The dashed line at 70 Mbar indicates the highest pressure 
explored in previous work. 



IA/mmm symmetry that appear at high pressure. More 
stable than C2/m above 135 Mbar, the P2\jm phase 
is also a monoclinic, layered structure that resembles 
the C2/m phase except that O-H-0 bonds are divided 
into two classes: one is shortened and straight, while the 
other is squeezed and distorted in different directions. At 
330 Mbar, P2\jm is replaced by a tetragonal IA/mmm 
phase, which has the same structure as the L'2 phase of 
ThH2 [19]. lA/mmm remains the most stable structure 
up to 500 Mbar. Like the C2/m structure, an analysis of 
the electronic band structure shows that the P2\/m and 
IA/mmm phases are metallic. Structural parameters of 
these two phases are listed in Supplemental Table I. 

In addition to our initial studies on water ice, we per- 
formed AIRSS on hydrogen (results see Fig. |2j. Below 50 
Mbar we find the same progression of ground state struc- 
tures, Cmca — > C2/c —> IAi/amd (or Fddd) — > Cracra, 
as has been reported previously [20H22] . When pressure 
exceeds 54 Mbar, we find that the 12-atom per unit cell 
Cmcm structure of Ref. [20 , hereafter Cmcm-12, be- 
comes less stable than a new Cmcm-A structure, in which 
three groups of bonds (different in length) exist without 



forming in-plane H3 clusters [20]. As pressure increases 
to 87 Mbar, the hep structure (symmetry PG^/mmc) be- 
comes stable. Different from the P63 / mmc structure pre- 
dicted in Ref. [21] (c/a « 2), our hep lattice exhibits a 
smaller c/a ratio of 1.58. The hep structure remains sta- 
ble until being replaced by the fee lattice at 180 Mbar. 

Having computed the ground state structures of H 
and H 2 throughout a wide pressure range, we now ex- 
plore novel hydrogen-oxygen stoichiometrics, including 
hydrogen-enriched structures H3O, H4O, H^O and HgO 
and oxygen-rich structures HO, HO2 and HO4. Between 
200 and 3000 random structures with up to 8 f.u. at 
pressures of 20-500 Mbar were generated, depending on 
the specific computational cost and searching efficiency. 

Figure [3] compares the stability of structures with dif- 
ferent stoichiometrics using a convex hull diagram [23J. 
This diagram shows the per-atom enthalpy of the com- 
pound H m O n , as a function of fractional concentration 
of oxygen, x = n/(m + n), relative to a combination of 
the two phases at the ends. A particular H m O n phase 
is stable relative to a mixture of two phases of different 
concentrations if its enthalpy lies below the line that in- 
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FIG. 2. (color online) Enthalpy versus pressure of different 
structures of hydrogen relative to IAi/amd [21] in the left 
panel and the 4-atom Cmcm phase on the right. The refer- 
ence for each phase is given in the legend. In the right panel, 
structures of the 12-atom Cmcm of Ref. [20] and our new 
4-atom Cmcm are compared. The dashed line at 50 Mbar 
indicates the highest pressure explored in previous work. 



terpolates between those two phases. In Fig. |3ja) we plot 
a convex hull diagram showing the enthalpies of possible 
H-0 phases at a pressure of 100 Mbar. The H4O struc- 
ture emerges as the most stable on the hydrogen-enriched 
side of the diagram, implying that any H-0 complex with 
an H:0 concentration greater than 2:1 will preferentially 
form an H4O+H or an H4O+H2O mixture, in contrast 
to the H2O+H2 mixture we expect at ambient pressure. 
On the oxygen-rich side of the graph, the HO, HO2 and 
HO4 phases are unstable relative to H 2 0+0 mixtures. 

In Fig. |3jb) we focus on convex hull for the hydrogen- 
enriched compounds including the H3O and H4O phases 
after having switched the right end from O to H2O. Besed 
on our calculations at 10 and 20 Mbar, we predict H 4 
to become stable relative to a mixture of H and H2O at 
14 Mbar. Such a pressure is predicted to for the core- 
mantle boundary (CMB) of Saturn and far below that 
of Jupiter. For pressures above 140 Mbar, we find H3O 
becoming stable relative to H2O+H4O mixtures. These 
results imply that at giant-planet CMB pressures, in the 
absence of temperature effects, excess hydrogen present 
is incorporated into existing water ice structures to form 
the novel hydrogen-rich compound, H4O. On the other 
hand, no oxygen-rich compounds will form, which is in 
agreement with recent calculations at high temperatures 

The H4O phase that we find here is a layered or- 
thorhombic structure with Prima symmetry. Each layer 
is characterized by -H-O-H-O- chains, of which each O 
is associated with one H 3 unit shown in Fig. [4j Our 
band structure analysis showed that the H4O structure 
becomes metallic at 80 Mbar. For H3O, we find a bee 
(Im3m) structure above 110 Mbar, a layered monoclinic 
C2/m structure (see Supplemental Fig. 1) at 20-110 
Mbar, and a slightly less stable orthorhombic Cmmm 
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FIG. 3. (color online) (a) Enthalpy of different H-0 com- 
pounds at 100 Mbar, relative to a mixture of pure H and 
O. (b) Enthalpy of hydrogen-rich H-0 compounds at several 
pressures, relative to a mixture of H and H2O. 




FIG. 4. (color online) Structures of H^O-Pnma (left) and 
H3O-6CC (right). Red and blue spheres denote O and H atoms, 
respectively. 



structure (see Supplemental Fig. 2) between 20 and 100 
Mbar, whose enthalpy is very close (2-36 meV/f.u.) to 
that of the C2/m structure. 

To verify the applicability of the PAW pseudopoten- 
tial under extreme pressures, we performed all-electron 
LAPW calculations using the ELK code [25 . The LAPW 
calculations for hydrogen at 80-180 Mbar reconfirm the 
Cmcm — >• hep transition at 87 Mbar and increase the 
enthalpy of the fee structure by 5 meV, precluding fee 
to become stable until a pressure of about 200 Mbar is 
reached. We also recalculated the relative enthalpy of 
two systems (H+H2O versus H3O, at 100 Mbar), and 
found that the difference between the two methods is 
only 4 meV/f.u.. We further replaced the VASP-PAW 
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pseudopotential of oxygen with an all-electron PAW po- 
tential (cutoff radius equals to 0.95 Bohr), re-optimized 
two H2O structures, IA/mmm and C2/m, and compared 
their enthalpy. The enthalpy difference, between the 
IA/mmm and the C2/m structures, increases by 16 and 
70 meV/f.u. at 180 and 400 Mbar, respectively. All 
these corrections are small and confirm our prediction 
of two new ice phases, although their transition pres- 
sures increase slightly (by < 4% for P2i/m, < 18% for 
lA/mmm). 

We also estimated the influence of zero-point motion 
(ZPM) by phonon calculations with 2x2x2 supercells us- 
ing the finite displacement method [26]. Including zero- 
point energy (see Supplemental Table II), the convex hull 
diagram in Fig. [3^b) shifts in the positive direction by 
12 meV/aotm for H4O and 3 meV/aotm for H3O at 20 
Mbar, and by 74 meV/aotm for H4O and 89 meV/aotm 
for H 3 at 180 Mbar. By applying these changes to the 
linearly interpolated enthalpy, we estimate the ZPM in- 
creases the pressure at which H4O and H3O become more 
stable than their corresponding H+H2O mixtures by, at 
most, 1.3 Mbar. Relative to the H2O+H4O mixture, H3O 
becomes less stable below 180 Mbar. We note that, based 
on the findings in Refs. [6j[8], the ZPM-induced changes 
in the enthalpy differences between H2O structures are 
small and do not change the sequence of phase transition, 
although the transition pressures will be slightly affected. 
Further, the difference in the zero-point energy of differ- 
ent hydrogen structures is small (< 80 meV/atom when 
below 50 Mbar), as indicated in Ref. [21 , and will not 
change our conclusions on the stability of H4O. 

In this work we have replicated and extended the phase 
diagram of solid hydrogen to 180 Mbar and that of wa- 
ter ice to 500 Mbar. We also predict novel, hydrogen- 
enriched meta-ice structures to become more stable than 
hydrogen- water mixtures at high pressure. Our results 
imply that water ice under conditions in excess of 14 
Mbar at low temperature will absorb excess hydrogen 
from its environment to form a hydrogen-rich H4O phase. 
The 14 Mbar pressure in question is comparable to 
the CMB of Saturn and far below the core pressure 
of Jupiter, and may potentially also be reached inside 
super-Neptune ice giants. This result stands in contrast 
to the tendency of methane to preferentially expel hy- 
drogen from its own structure under pressure to form 
hydrocarbons and eventually diamond [27] . and in ice 
mixtures methane could potentially provide the excess 
hydrogen to incorporate into H4O. These results under- 
line the fact that chemistry at high pressure may deliver 
substantially counterintuitive results, and that consid- 
eration of structures likely to form at giant-planet core 



conditions requires looking beyond traditional ambient- 
pressure chemistry to explore unfamiliar stoichiometrics 
and combinations of elements. 
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